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Abstract: We present a new systematic method to evaluate one- loop tensor inte- 
grals in conventional ultraviolet cutoff regularization. By deriving a new recursive 
relation that describes the momentum translation variance of ultraviolet integrals, 
we implement this relation in the Passarino-Veltman reduction method. With this 
method, we recalculated the Higgs boson decay into two photons process at one-loop 
level in the Standard Model. We reanalyze this process carefully and clarify some 
issues arisen recently in cutoff regularization. 
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1. Introduction 

Recently, the ATLAS |l| and CMS[[| collaborations have renewed efforts to search 
for the Higgs boson at the CERN LHC with data integrated up to 0(fh~ 1 ). The 
excluded mass region for the Standard Model(SM) Higgs boson has been extended 
to most of the region between 145 and 466 GeV. In the low mass region of the Higgs 
boson, the two-photon mode of Higgs decay plays a crucial role in experimental 
studies. 

R. Gastmans et al. recently recalculated H — > 77 via W-boson loop |f| f§], 
which yielded a result in contradiction with the old ones in the literature [|, |], [7], §]. 
Their computation was carried out in four-momentum cutoff regularization rather 
than dimensional regularization (DREG). To reduce the number of Feynman dia- 
grams, Gastmans et al. chose unitary gauge. In their treatment, the new result, 
which satisfies the decoupling theorem^, was favored by the authors. Later, sev- 
eral authors |L0], [LI], 0, [L3| have pointed out that the old results are still correct 



and decoupling theorem is violated by Hcj) + (j)~ in this case. However, we are still 
unsatisfactory with the explanations about the problems with the calculations of 
R. Gastmans et al., since there have never been doubts about the correctness of their 
algebra. In order to clarify this problem, we develop a new method to do one-loop 
calculation in cutoff regularization. 

Although DREG has proven its superiority and achieved the most widely usage 
in phenomenological applications, cutoff regularization, the oldest regularization, still 
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has some advantages compared with DREG theoretically. For instance, in DREG, 
one is unable to obtain the correct divergent terms higher than logarithmic diver- 
gences, which means that quadratic divergent terms of SM Higgs self-energy diagrams 
disappeared in DREG. Pauli-Villars regularization is flawed because it violates chi- 
ral symmetry, while the symmetry is preserved in cutoff regularization. Moreover, 
from Wilsonian effective field theory viewpoint, cutoff regularization scheme is also a 
more intuitive and straightforward scheme. Therefore, the introduction of an explicit 
cutoff is sometimes advantageous. 

However, there are still many drawbacks in this four-momentum regularization 
that should be mentioned. Considering the truncation in momentum modes, this reg- 
ularization is flawed because it violates gauge invariance and translation invariance 
regarding the loop momentum. The latter condition signifies that the results may 
ambiguously depend on the manner how the propagators are written. Hence, in the 
present paper a new recursive relation for loop momentum translation is derived first. 
Then the Passarino-Veltman reduction method 1 [1141 Il5l is modified to reduce the ten- 



sor integrals in this regularization. One can follow Dyson's prescription |L6| [L7| to 
obtain a gauge invariant result, just as shown in our calculations of H — > 77. 

As an example, we reconsider the process H — > 77 in this four-dimensional 
momentum cutoff regularization with our proposed approach. Given that the process 
of H 77 is free from infrared and mass singularities, only the ultraviolet cutoff is 
considered here. Readers who need to handle infrared or mass singularities should 



turn to the mass regularization scheme demonstrated in the literature e.g.|18 

The present paper is organized as follows. In Section ^ a new recursive relation 
for the loop momentum translation in cutoff regularization is demonstrated. Then, 
it is implemented in the Passarino-Veltman reduction schemes in Section [3]. With 
this approach, calculations and analysis of H — > 77 are performed in Section f|. Our 
conclusion is present in Section |[ In Appendex|A|, the expressions for used in 

Section |2| is derived. Finally, some scalar integrals can be found in Appendix |B|. 

2. A New Recursive Relation 

In this section, we will show how to calculate 



lll...fls 



'-{k - b) 2 + a 2 ) n J (-k 2 + a 2 ) n 

This integration have a superficial divergence degree A = s + 4 — In. 



1 Note that, integration by parts (IBP) reduction methods are not valid in this case due to 
nonvanishing surface terms. 
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A negative A evidently simplifies the calculation, because the limits of the inte- 
grals in Eq. (|2.1| ) can be set to infinity and the translation shift k — )■ k + b does not 
change these limits. Therefore, 2r Ms completely vanishes when A < 0. However, re- 
sults may vary when the integrals in Eq.( [2.1| ) are ultraviolet divergent because there 
is an artificial four-momentum cutoff scale A in these integrals. These conditions are 
then considered following. 

1^ ,, can be rewritten as 

LL\ ...Us 



tA (h n 2\ _ ( f j4. Kl ■■■ k ^s _ f ,4, k^... k^ \ 

[J dk ^ {k - b) 2 + a 2 )n J d \_ k 2 + a 2 )n ) 

with the remainder f rem {b) = (k + 6) w . . . (k + b)^ — k^ . . . k^ s . Using the identity 



1 1 - f\ x < B ~ A ) (23 ) 
A n B" J {xA+{l-x)BY+^ K ] 



one arrived at 



<U(M 2 ) = I d'k I dx 



i n {-2b- k + b 2 )k^...k^ a 



{-(k - c) 2 + d 2 ) n+1 



-k 2 + a , 

where c = x b, d 2 = a 2 — 6 2 x(l — x). The integral momentum to k + c in the first 
integral of the previous equation is shifted, so that 



Jo Jo 



,4,„ I A n(-2b-k + b 2 -2b-c)k, 1 ...k, s 



+ d A k dx 



(-k 2 + d 2 ) n + l 

+ /d , fc /" d ^-*+'' 2 - 2 ''-)/" m w 



-k 2 + d 2 ) n + l 



- / d*k 



-k 2 + a 2 ) n 
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-1" / /'dx/^Ccd 3 



o 



-2;, / d'A" / dx k ^K,---Ks 







-A; 2 + d 2 )^ 1 



-2n I d*k I dx 1 



K -k 2 + d 2 ) n + 1 

df rem {c) 

d 4 A; / dx-—^—-. (2.5) 
J (-k 2 + d 2 ) n v ' 

In the sixth line of Eq.( |2.5| ), a spurious part that is proportional to (1 — 2x) in 
the numerator of the integrand is removed. What's more, integration by parts is 
performed at the end of Eq. (|2.5|) . 

Aside from the terms expressed in J A_1 and I A ~ 2 , Eq. fl2.5|) can be simplified 
further using the formulae ^...^(a 2 ) = / d 4 fc -A^j^g given in appendix After 
expanding the terms proportional to a;- 7 in the integrands and implementing the 
expressions for J , a lot of terms are canceled. Thus, the final result is 



/A ( 6 a 2) = _2n 6 W) dx/^- 1 u (c, d 2 ) +n6 2 / dxI A - 2 (c, d 2 ) 

fll.../l s \ J / I IMJUL.-Hs \ ' / I fj,l...fl s \ ' y 

JO Jo 

t=max(0,4— 2n), even ^2 / 

(2.6) 

where the notation {5 ,s+ *~ A & A ~*} ( u 1 ... jUs defined in appendix [A], n = s+4 j~ A , |_j/J is a 
Gaussian function ( the greatest integer that is not larger than y), and the function 
is 



a na,ni,no>0 / -|^ 2 +n n fA f) ;| 

V 7 f-f ^ f 2n 2 + ni + A-tn 2 !n 1 !n ! V 7 V 7 

t=0 ri2+rai+no=( 

(2.7) 

We should make some remarks about above equation before going forward. In the 
recursive relation Eq. (|2.6| ), there are two integrals left. However, since all of the 
I A ^ (6, a 2 ) are only polynomials of b and a 2 ,i.e.,they are only polynomials of integral 
variable x, the explicit expressions for this recursive relation can be easily obtained 
with the help of computers. Especially, the explicit expressions for J A „ a (A = 
0,1,2,3) are 
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r( 2 , 

ll...Jb,a 2 ) = (n b 2 {g°}, 1 ..., s + (4n + 2)6(s-2){g s - 2 b 2 }, 1 ..., s ) 

i?r 2 (-2) 1 - n , s + 2 

— ; , with n = , 

r(n + 2) ' 2 ' 

,3 r . 2x ^^ 2 (-2) 2 -(3n 2 + 6n + 2) 

15 j pi. ..[is r(n + 2) 



A - n a i -b 2 , with n = . (2.8) 

n + 2 J ' 2 v 7 

3. Modified Passarino-Veltman Reduction Schemes 

It is known that the one-loop tensor integrals can be reduced to a linear combination 
of up to four-point scalar integrals [Q. In this section, a generic one-loop integral 



T, 



N 



d A k- 



k k 
^Hx ■ ■ ■ ^fis 



(3.1) 



D A . . . Av_i 

with propagators = (k + Pi) 2 — m 2 + ie and po = is considered. As mentioned 
in the previous sections, the integrals ^ s may not be translation invariant be- 
cause of the finite integral limits when they are ultraviolet divergent. Therefore, the 
expressions for 



AL 



N 



ki. 



(k-Pi, 



D 1 ...D N J D 1 ...D N 
{k+pi) 2 — mf+ie but pi ^ and Di = (k+pi 



(3.2) 



pi) 2 -m 2 + ie 



where propagators D % 
should be calculated. After the conventional Feynman parameterization, AL^ ^ s 
can be reexpressed as 



N 



simplex j_i 



+4,-2N, 
Ms \ 



Pi, a 2 ) 



N 



simplex 



s-i(pS-iji+4-2N, i( h ~2 



(M )} w . 



i=0 



, (3.3) 
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where b = - UiPi + Pi, a 2 = - J2?=i u i(Pi ~ m i) + Z*j=i ^i u jPi ' Pj, and the 
notation {pl~ l P +4 ~ ,% (b, a 2 )}^...^ is defined in appendix [A| with the divergence 
degree of p +4 ~ 2N ' 1 defined in Section |] is i + 4 — 2N. The polynomial dependence 
of b, a 2 in I A (b, a 2 )^...^ makes the simplex integration in Eq. (|3.3|) straightforward 
using 



./simplex i=1 i=1 1 {l^i=l r i) 

Next, several notations similar to that given in ref.||15|| are reintroduced here 



no,ni,...,n.jv>0 

ALjL* = {9 2no P? ■ ■ ^} W .., S A/>L,.„ : 

2n +™l+---+"jV=S 2n rijy 

rao,ni,...,?ijv_i>0 

— Li/ Pi ■ • -FAT-1 i/-ii---Ms J 0...0...(V-l)...(7V- 1)> 

2n +ni + ...+n JV _i=s 2n„ " TTT - ' 



^//.i • • • kfi s 



Tl..Jk) = J d*k 



N 

k k 



D ...D k ...D N 

k k 



D X ...D N 

no,ni,...,npf>0 

2n +n 1 +...+n N =s 2nQ npf 

rao,ni,...,?ijv_i>0 

^(0) = £ {9 2n °(P2-Pi) ni • • • (PN-Px)""- 1 }^ 

2no+ni + ...+ri]v-i=s 

^0. . .0... (N — 1) . . . (N — 1) 

(0), (3.5) 

2n n JV _ 1 

with D . . . D k . . . D N = D . . . D k _iD k+ i . . . D N and employing the caret " ~" to 
indicate the indices omitted. 

Thus in this cutoff regularization, Eq.(2.9) in ref.fjTSj should be replaced by 

k AT-1 

^0^1_^_l,i 2 „ +fc+1 ...i s (0) = (~) k ^2C l k 2j ^o)...,0,ii,...,ii, i2n+fc+i-l,-, 
2n fc i = ii,...,ij = l " ^ " 

+ AL o^^ J2n+fc+1 ..i 3 > Wfc+i> • • > !• (3-6) 

2n fc 



- 6 - 



When the determinant of Gram matrix 



/ 2pi • pi ... 2px ■ p N \ 



Z (N) 



(3.7) 



\ 2p N - pa ... 2p N ■ p N J 
for (N + l)-point functions is non-vanishing, the reduction can be continued as 



rpN 

1 00i 3 ...i s 



2(3 + s-N) 



N-l 



T: 



N 

ii...i s 



3=1 

3=1 

2V-1 / 



N 

iz...i s 



N-l 



' 00l2...ir— is 



, ii^O, (3.8) 



3=1 



r=2 



where some notations are defined in ref. 15 



r _ 2 2 i 2 



Urjfc = 



5S 



T, 



N-l 



Xk)8 ki2 ..J kia -T»:l(0)-f k Ti 



iN 

i 2 .. -i a ' 



(3.9) 



ki 2 ...i s J- (i 2 ) fe ...(j s ) fc v 

Otherwise, when the Gram determinant is zero, there is at least one non-vanishing 

- 7 {N 
J kl 



element Z9P in the adjoint matrix of 



L ki 



:-) 



k+l 



2pipi 

2pfe-ipi 
2pfc+iPi 

2p N Pi 



2p 1 p l _ 1 2pipi + i . . . 2p x pN 

2p fc _ip/_i 2p fc _ipi+i . . . 2p fc _ip7v 
2p fc+ ip;_i 2p k+1 pi +l . . . 2p k+1 p N 

2p N p l _ 1 2p N p l+1 ... 2pjvpjy 



(3.10) 



and one non-zero element Xq- in the adjoint matrix of the following matrix 
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X (N) 



(2m\ fi ■■■ fN \ 
fx 2ptpi . . . 2p x p N 



\ f N 2p N p x . . . 2p N p N J 
One-loop reduction can be applied using the following equations 



(3.11) 



rpN 

ii-.-is 



rpN 

1 00i 1 ...i s 



N-1 



^(N-l) in \ °ni\...i a * Z^ ™ r 00ii...i r . 

r=l 

\y(N-l) qs+2 

y^ki J mi...i s 



(N-1) 
Qj n=l 



2(6 + ,-iV + E: =1 4o)^f- 1) 



N-1 



' y^nl °nki 1 ...i s n kl °nni\ 



n=l 

N-1 



n,m=l 



r=l 



f nf mXi 1 ...% a 2 ^ ^ {fn^mi r fm^ni r ) 



m00ii...i r ...is 



N 

00ii...i r ...i s 



r=l 



->N 

' OOOOii... i r ... it ■■ -i s 



4 ^ 5 n i r 5 m i t T Q( 

r,t=l,r^t 

Some notations in Eq. (|3.12| ) should be recalled, i.e. 



(3.12) 



UN) 

'mm 



\i+j+k+l 



7(2) 

OS 



sgn(z — &)sgn(/ — j) 
2p x pi . . . 2p x pj-i 2pip j+1 . . . 2pipi_ 1 2p x pi + i 



2pi-iPi 
2pi+iPi 

2pk-iPi 
2pk+iPi 

%PnPi 
N>2 



2pi_xPj-i 2p^ip j+ i 
2p i+1 p j ^ 1 2p i+ ip j+ i 

2p k -ipj-i 2p k _ip j+1 
2p k+1 p j _i 2p k+ ip j+1 



ZpNPj-i ZPnPj 



+1 



T {i2)l..{is)S k ^ ki * 



2pi-ip«-i 2p^ipi + i 
2p i+1 p^i 2p i+ ipi +1 

2p k _YPi-i 2p k _ipi + i 
^Pk+iPi-i 2p k+1 p l+1 

2p N p l _ 1 2p N p l+1 



r rpN-1 



(0). 



2pxp N 

2p i _ 1 p N 
2p i+1 p N 

^Pk-lPN 
Zpk+lPN 

2p N p N 



(3.13) 
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For det(Z( 7V - 1 )) = 0,det(X( 7V ~ 1 )) = and all X { 
X^-~ l 7^ 0,following equations 



{N-1) 
Ofc 



but Z ( k f' 1] ^ and 



T, 



N 

. . . 1 . . . li\...i m 



r n 



T; 



N 

i\...i a 



2(n+l)^f- 1} 



O y{N-l) T N 

/ j hU O...Ol...li 1 ...i j ...i rl 

7- n+1 



^0 1 . . . lii...i m 

r — 2 n + 1 



' ^1 j • • • ) 7^ 0, /, 



N-1 

_|_ giN-l) gr+n+r 

/ j kj j0 ■ ■ 

3=1 

1 ^^r i H 2 (5+^-^« i ,., s -^i(o) 



(N-1) 



N-1 



n=l 

AT-1 / 



m,n=l 



r=l 



(3.14) 



can be used. Other details in the derivation of these equations can be found in [|15 



4. Higgs Decay into Two Photons 

In this section, one-loop reduction as illustrated in the previous section is applied to 
the process H — > 77. 

In unitary gauge, the three diagrams via the W-boson loop that contribute to 
this process with a specific loop momentum configuration are shown in Fig. ([!]). A di- 
rect calculation of amplitude (dropping the polarization vectors of external photons) 
yields 



" unitary 



with 



3e 3 mw 



[2k%k" (in 2 - (m 2 H - 2m^)C (0, 0, m 2 H , m 2 w , m 2 w , m 2 w )) 



-h ■ k 2 g^ u (in 2 - 2(m 2 H - 2m^)C (0, 0, m 2 H , m^, m^, m^))] 
3ie 3 mw 



k^kl 2m 2 H + 4(m 2 H - 2m 2 w )f( 



Am 2 w 



+h ■ k 2 g^ m 2 H + A{m 2 H - 2m 2 w )f ( 



m 



1-1 



Am 2 w - 



(4.1) 



f(aO = 



arcsin(y / x) 2 
lnf 1+ Vl ~ x ~l ) - 



1 2 



- in 



,X<1 
, X > 1 



(4.2) 
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and the scalar integral Cq is given in appendix [B| However, gauge invariance is 
spoiled in this four-momentum cutoff regularization. Therefore, a term should be 
subtracted from the above expressions to recover gauge invariance. In this gauge, a 
requirement of M.^ v (ki = k 2 = 0) = should be made. However, 



J unitary 



(fcx = k 2 = 0) 



-3ie 3 mw 

167T 2 S„, 



<r ± o. 



(4.3) 



Following Dyson's prescription [|16|, |I7|, gauge invariance is recovered after making 
subtraction from Eq.([4.3[), and the final result is 



with r = tHt- following the notations of refs.[§, f|]. Eq. ( |4.4|) is the same as those 
in refs. |3], f|] up to a factor of — 2i from the symmetry factor of loops and different 
conventions of Feynman rules. However, in this gauge, there are high degrees of ul- 
traviolet divergence in each diagram. The expressions for amplitude may be different 
under different choices of loop momentum. One may suspect that the discrepancy 
between Eq. ( |4.4j) and the result given in DREG 



jp3 

Kbbg = o - - sTh ■ **) (2 + Sr- 1 + 3(2r- 1 - r~ 2 )f (r)) (4.5) 

107T //i^^/Oy; 

is originated from the bad loop momentum choices in Eq.( (4.4| ). However, in our 
calculation we find that the terms 



ie 3 



AM^{p) = -jr^-z [{Kp v - V»K) (-3A 2 - 2m 2 H - Qm 2 w + {k x + k 2 ) ■ p - p 2 ) 



+2g^(h -k 2 )-p (-3A 2 - 2m 2 H + 3m 2 w + (h + k 2 ) ■ p - p 2 )] (4.6) 



should be added to Eq.( [4.4| ) if loop momentum k is shifted to k + p. From the 
symmetric consideration of k±,k 2 , a proper choice of p is fcl + fc2 which is the same as 
that presented in refs.g, g. Since AM^ U {^^) = in Eq.([0|, the result in Eq. Q 
remains unchanged. From Eq. (|4.6j ), it seems hopeless that the difference between 
Eq. (|4.4|) and Eq.(4J5) can be eliminated through shifting the integral momentum k. 



In 't Hooft-Feynman gauge (£ = 1), the amplitude with one-loop diagrams shown 
in Fig. (|) is 
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+Qm 2 v {m 2 H - 2m 2 v )C {0, 0, m 2 H , m 2 w , m 2 w , m^)] + hi ■ k 2 g^ u [i-n 2 {m 2 H + §m 2 w ) 
-12m 2 w (m| - 2m 2 w ) C o (0, 0, m 2 H , m 2 w , m 2 w , m 2 w )] } 



ie 3 



[-2k\k% (mjj + 6m 2 H m 2 w + 12m 2 w [m\ - 2m 2 w ) f (r)) 



16n 2 m 2 H mwSv 

+h ■ hgT {m A H + 6m 2 H m 2 w + 2Am 2 w {m 2 H - 2m 2 w )f(r))] . (4.7) 

Following a similar procedure to obtain a gauge invariant result, the amplitude at 
k\ = k 2 = is calculated as 



which is the same as Eq.( f4.4|) . However, the gauge invariant amplitude is non- 
vanishing at k\ = k 2 = because of the contributions of diagrams (g) and (h) in 
Fig.(0) in this gauge. These contributions are 



M^ {gM (k 1 = k2 = 0) = ^ I -r. (4.9) 

Therefore, the subtracted terms should be Mf^L-^ki = k 2 = 0) — M^Li i g h \ (k% = 
k 2 = 0) instead of M.^Li{k\ = k 2 = 0). The final result is 



M t=i = -77-2 - 9» u h ■ k 2 ) (2 + Sr- 1 + 3(2r- 1 - r~ 2 )f (r)) 

= ACrbo- ( 4 -!0) 

The term generated by the contributions of the Goldstone triangle diagrams (d, e) 
in Fig.(0) spoils the decoupling theorem, as pointed out by Shifman et al. recently 



1C]. Given that there are only logarithmic divergences under this covariant gauge, 



the result in Eq.( [4.10D is unique with a different loop momentum chosen. To best 



of our knowledge, it is the first derivation in the 't Hooft-Feynman gauge in cutoff 
regularization. 

It seems there are some problems with unitary gauge in this cutoff regularization. 
Given that the top quark loop ( Fig.@) does not suffer from any ambiguities in the 
gauge or loop momentum choices, the diagrammatic expressions are expected to be 
the same in DREG and in this cutoff regularization. These conditions have been 
verified following the same procedures. The result is as follows 
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AC 



ie 3 N r 



187r 2 m w s w 



(k$k» - h ■ k 29 n (x' 1 + (X- 1 - X- 2 )f(x)) , (4.H) 



4m? 



where x 

The authors of refs. flll], |12| have also calculated this Higgs decay process in Pauli- 
Villars regularization and dimensional regularization respectively, and obtained the 
same result as the old ones|5|, ^ 0, |j. Their statement about this issue is that the 
integral (in Euclidean) 



is vanishing in cutoff regularization, while it is nonzero in DREG, which is also 
pointed out by R.Gastmans |3|, Q. They also argued that contained the difference 
of two logarithmic divergencies and should be regulated. Therefore, the integral that 
violats electromagnetic gauge invariance may suffer from some ambiguities. Actually, 



this issue was first discussed by R. Jackiwpifl in a more general case. However, we 
think that the vanishing of in 4 dimensions is just a result of the fact that the 
integral intervals are symmetric about the origin even when there is a cutoff A, and 
a replacement of k^k v — > 9p,u . in the integrand is also proper. 

Moreover, very recently R.Jackiw also pointed out that by combining the two 
terms in the integrand of one can avoid infinities but the difference of the integrals 
in these two regularization schemes is still the same, thus both evaluations are math- 



ematically defensible |22|]. So, what are the physical reasons for these ambiguities? 
In the following, we will try to clarify this issue. 

Considering that the diagrammatic expressions in unitary gauge are not well- 
defined in the four-momentum cutoff regularization, how to recover the correct result 
under this condition may be still an open question. We investigate the Lagrangian 



for the Standard Model in unitary gauge, similar to the treatments in ref.[|19|]. The 
covariant terms for scalars are 



^scalar = (D^^D^) - V($), 

with V($) ee -fi 2 &$ + A($ t $) 2 ,D At $ = (d„ - ^gr'W; - |g'B^ $. (4.13) 

The scalar doublet produces the vacuum expectation value through the Higgs mech- 
anism as 



- 12 - 



\ f^ 1/2 



V2 

Therefore, the scalar fields can be redefined as 



<*>«= I v ],v=|y) ■ (4-M) 



There are terms like 



m^W+W"'' 1 + im w {W~d^(f) + - Wj<9^~) 



rn 



w 



m w 



W 



mw 



after expanding the Lagrangian given in Eq.flOSD, where W± = ^ M . By follow- 
ing the prescription in ref . [119 , the W-boson fields in unitary gauge can be redefined 

as 



(4.16) 



* * m w 
W-=W-- —8^-. (4.17) 

In this gauge there are no kinetic term c^+c^ - and mass term for the Goldstone 
+ , 4>~ because of the cancelation between the last term in Eq. (|4.16 ) and the original 



kinetic term of the W-boson's Goldstone in £ sca i ar . However, terms such as h(p + (p 
still exist in the original Lagrangian. In DREG, <p + = <fi~ = can be set safely, similar 



to a previous work by Grosse-Knetter | 20fl , because all the momentum modes can be 
included in this regularization 2 . Hence, the conventional Lagrangian in unitary gauge 
only with physical fields is obtained. However, the results are in contrast to those of 
the four-momentum cutoff regularization, because an artificial scale A is introduced 
in the Lagrangian. The absence of a kinetic term for (fi + , does not mean that these 
Goldstone fields are vanishing intuitively, but because the theory does not provide 
any information above A in this regularization. If the mass of <f> + , <fi~ is assumed to 
be 0(A), there are still finite contributions from the Goldstone triangle diagrams 
when A — > oo. From this viewpoint, the cutoff regularization in unitary gauge is 



2 Note that the limits of loop integrals are taken to be infinity 
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( a ) 0) (c) 

Figure 1: Feynman diagrams via W-boson loop in unitary gauge for H — > 77. 

problematic. The violation of the property in gauge invariance can also be attributed 
to the absence of large momentum modes. Therefore, the old results in the literature 
for H — > 77 are still valid. 

In fact, dimensional and Pauli-Villars regular ization schemes are free of missing 
large momentum modes, and can maintain gauge invariance. Therefore, these results 
are correct also in unitary gauge. From the evaluations of R.Jackiw, the integral 
can be dealt without any infinities, and the only difference is from surface terms 
(i.e., large momentum region), which also verifies our conclusion. 



5. Summary 

A method for systematical evaluations of one-loop tensor integrals in cutoff regular- 
ization is proposed by deriving a new recursive relation Eq. ( j2.6|) and implementing 
it in the Passarino-Veltman reduction method. The result has been expressed in 
a form that can be directly translated into computer codes. Similar to the meth- 
ods presented in ref.[[H|, our results are also numerical stable for up to four-point 
integrals. Surely, our method can be extended to deal with high-point integrals 
straightforwardly. 

With this approach, we have calculated the amplitudes for Higgs decay into two 
photons via the W-boson loop and the top-quark loop. The correctness of the method 
has been confirmed by evaluating these processes and checking other programs, and 
it is certainly useful in both theoretical and phenomenological aspects. Moreover, 
we also reanalyze the Higgs decay process and make our efforts to find the physical 
reasons for some puzzles appeared in the calculations of this process. 
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(9) (h) (») 



Figure 2: Some representative Feynman diagrams via W-boson loop in 't Hooft-Feynman 
gauge for H — > 77. 




(a) (b) 



Figure 3: Feynman diagrams via top-quark loop for H — > 77. 
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A. Derivation of Expressions for J^ lfl2 fls 

In this appendix, the general formulae for </^ 2 ... Ms = / d 4 fc ^p^w used in section 
H are derived first. Obviously, Ju llMl ...u a is vanishing until s is even. Therefore, s 
should be an even and non-negative integer and N should be positive in the following 
context. 

A notation (similar to but a little different from that in Ref. ||15||) is introduced 
first in order to write down the tensor decomposition in a concise way. We use 
curly braces to denote symmetrization with respect to Lorentz indices, where all non- 
equivalent permutations of the Lorentz indices on metric tensor g and momenta p are 
implicitly understood. A generic notation {g^p™ 1 . . •Pk k }m...fH with t = Y^i=i n i + 
2n means a sum that the 2n of Lorentz indices fix, . . . ,m are distributed to n 
metric tensors g while ni of them are distributed to ni momenta pi with equal 
weights. For instance, 



p,upa — 9pv9pa 9pp9v(J 9p,a9vpi 

{9 2 p 1 }pup = 9pvV P + 9v P Vp + 9ppPv, 

{PlPD^P = P\pV\vP2p + P\pP\pP2v + PluPlpP2n- (A.1) 

The Lorentz covariance ensures us to make the following replacement 



kp, 1 kp l2 ...k fls > {g }p ll p 2 ... fla — 2^2% (A-2) 

in the integral , which can be proven by the induction of the integer s. After 

this replacement and subsequent Wick rotation, spherical coordinate system trans- 
formation and some trivial variable substitutions, one arrived 



-A 2 ■ - - - 

r(f + 2) J ^{K + a 



T» inn ^ 2 (-2)^ [ AfC K2 (L* 



where A was denoted as the ultraviolet cutoff scale. 
Eq. (|A.3|) can be solved directly when a 2 = 0, i.e. 



J pi...p s ~ 19 jpi...ps _|_ 2) A i^-j; 
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when superficial degree of ultraviolet divergence A = s — 2N + 4 > 0. In the case 
of A < 0, Eq.( A.3 ) encounters infrared divergence, which is not considered in this 
article. When a 2 ^ 0, result becomes a little more complicated than the previous 
one. However, this problem can be resolved after implementing the tricks of using 
integration by parts and following integral formulae 



/I / ™ ( _„\k n+l-k 

dx x n ln(x + a) = — — x n+1 ln(x + a) - V ^— '— — — 
n + 1 V f—f n + 1 — k 

-(-a) n+1 \n(a + x)) , n G N (A.5) 

into Eq.( |A.3| ). The explicit expressions for J^__ can also be obtained, i.e. a 2 ^ 
and superficial degree of divergence A = s — 2N + 4 > yields 



J, 



N 

Hi... Us 



{g 3 }^ (-2)- 



r(N) 



N-l 



T(N-k) 



£r(l + 3-^) 



rfi 



^ + D I ^ 

2 ^ L ) \ k =0 



—a 



2\k 



2A 



A-2fc 



A - 2k 



+ (-a 2 )2 ln( 



A 2 



(A.6) 



while a 2 ^ but A < returns 



JL, S = tfW-M. (-2)-f ^^a A (A.7) 
B. Some Scalar Integrals 

After the reduction of one-loop integrals using the modified Passarino-Veltman re- 
duction formulas given in the section |2|, every tensor integral can be expressed as a 
linear combination of up to four-point scalar integrals. In this appendix, the ana- 
lytical expressions for some scalar integrals are listed below. Some of them may be 
used in the one-loop calculations of the process Higgs decay to two photons. 

First of all, the conventions for the scalar integrals used in this article are fixed 
as follows 



For one-point functions, 



A (0) = -in 2 A 2 , 

( A 2 A 2 
A (ml) = m 2 m\ ln(— ) f ) , 



= r(ir^" (E U — A-,„r 2 - 2i + H^CM^) )b.2) 



Two-point functions can be easily verified as 



B (pl ml m\) = itt 2 ( ln(^) + 1 + ]T [ 7i ln(^-^) - ln( 7i - 1)] 



p 2 — m\ + m§ ± \J (p\ — m\ + ml) 2 — Ap>\m\ 

with 7 i 5 2 



2p 2 



A 2 

B (0,0,m 2 ) = ^ 2 ln(— ), 

m z 



fi (p 2 ,0,0) = m 2 I ln(Ar) + 1 



p2 



Bob 2 , 0, m 2 ) = in 2 ( ln(4) + 1 + In^ 



' m 2 p 2 m 2 



m„ ln(Ar) — m? ln(- 

D / n 2 2\ . 2 u Vm - y 1 

-dq(0, rn , m x ) = m — 



2 2 

mg — m| 



B (0, m 2 , m 2 ) = in 2 (ln(^r) - l) • (B.3) 



Moreover, two special finite three-point functions are 



4m 2 

C (0,0,p 2 , m 2 ,m 2 ,m 2 ) = ^ ( ln 2 f V / L ) - tt 2 



C (0,0,0,m 2 ,m 2 ,m 2 ) = -^. (B.4) 
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